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A BIHARMONIC RELAXATION METHOD FOR CALCULATING THERMAL STRESS IN COOLED 

IRREGULAR CYLINDERS^ 

By Ahthce G. Holms 


SUMMARY 

A numerical method vxls developed Jor calculating thermal 
in irregular cylinders cooled by one or more internal 
passages. The use oj relaxation methods and elementary meth- 
otis of finite differences was found to give approximations to the 
correct values when compared with previously known solutions 
for concentric circular cylinders possessing symmetrical and 
asymmetrical temperature distributions. 

INTRODUCTION 

USE OF COOLED IRREGIXAB CTUNDERS 

The evolution of aircraft propulsion systems has led to the 
frequent employment of cooled structures. The concentric 
hollow cylinder is a fanuliar example, although in some cases 
irregular cylinders such as cooled turbine blades with several 
internal passages (reference 1) are under consideration. A 
method of calculating stresses in thin-waUed turbine blades 
of the air-cooled type is presented in reference 2, but the 
general problem of calculating thermal stresses in long, 
hollow, thick-waUed irregular cylinders has not been solved. 

PREHOUS WORK 0.\ THERMAL STRESSES IN HOLLOW CTUNDERS 

Several methods of calculating thermal stresses for various 
special distributions of temperature in long hollow cylinders 
of particular shape have been developed. Some theoretical 
aspects of the general problem have also been discussed. In 
reference 3 the problem is regarded as an ordinary stress 
problem with given body and surface forces replacing the 
effects due to temperature distribution, whereas in reference 
4 the equilibrium and boundary conditions of the theory of 
elasticity are used without modifications to exhibit the tem- 
peratiue effects as body and surface forces. The method of 
reference 4 was applied to several special problems that had 
already been solved as well as to problems of composite 
bodies and eccentric circular cylinders. No application 
of analytical or numerical methods of calculating thermal 
stresses in cooled irregular multiply connected C3’linders has 
been published. 

SCOPE OF PRESENT I.NVESTIGATIOX 

An investigation was conducted at the NACA Lewis lab- 
oratory during 1949-50 to calculate thermal stresses in cooled 
irregular multiplj’ connected cj’linders. The problem of ther- 
mal stresses in irregular cj'linders is formulated in a manner 


that permits solution bj' the use of finite-difference methods. 
The contour integrals of reference 5 expressing the sin^e- 
valued character of the displacements for arbitrary' circuits 
around the internal boundaries are written in forms suitable 
for numerical methods of differentiation and int^ration. 
The boundary conditions based on the assumption of force- 
free boundaries and single-valued displacements are formu- 
lated in terms of derivatives of the stress function as sug- 
gested for uniform-temperature problems in reference 5, and 
stress functions are set up in a manner that is an extension 
to the thermal-stress problem of the work of reference 6 on 
doubly' connected domains at uniform temperatures. The 
relaxation techniques of reference 7 are used in solving the 
finite-difference problem of determining the stress functions. 
Details of the method are illustrated by examples. The 
method is applied to a symmetrically heated, hollow circular 
cylinder and also to a hollow circular cylinder with asym- 
metrical heating to show that the relaxation technique gives 
approximations to exact answers obtained by direct math- 
ematical methods. Comparison of stresses calculated by 
the relaxation techniques with those determined by exact 
methods is also made to compare the relative importance of 
several sources of error. The calculations for the concentric 
cylinder are described in detail sufficient to permit the method 
to be applied to more insular cylinders. 

SYMBOLS 

The following syrmbols are used in this report: 
an,«fl, 0(3 constants of integration 

E modulus of elasticity in tension and compression 

l,m direction cosines of normal drawn outward from 

region bounded by plane curve 
distance in zy-plane normal to plane curve 
residual in relaxation calculation 
polar coordinates 

arc length of plane curve in ry-plane 
temperature above initial stress-free state 
components of displacement 
rectangular coordinates 
coefficient of linear thermal expansion 
shearing strain components in rectangular coordi- 
nates 

unit elongations (straios) in x-, y-, and z-direc- 
tions, respectively' 

Poisson’s ratio 


n 

Q 

r,d 

s 

T 

u,c,w 

x,y,z 

a 


‘ Supersedes NAC.A TN 2iS4. "A Bihartnonfc Relaiatloa Method for Calculating Thermal Stress in CooM Irregular Cylinders" by Arthur O. Holms, 1851. 
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<rx,o-„,<rt aoraial components of stress parallel to x-, y-, 
and 2 -axes, respectively 

Txg sheaiing-stress component in rectangular coordi- 

nates 

<j) Airy’s stress function 

« component of rotation about s-axis 


where is Airy’s stress function defined by 


* by* 


L 

‘"■'"da:* ( 


(2) 


V* 


V* 


harmonic opei'ator 




biharmonic operator ^+2 


5* d* 


THEORY 

ASSUMPTIONS 


bV 

dxby J 

As shown in appendix B, the axial stresses can be calculated 
from the equation 

fft=E {ax-}-by-]rc)'\-v((rx-\-ff^ — aET (3) 


The following conditions are assumed to exist: 

1. Steady-state heat flow exists with temperatures uni- 
form along any line initially parallel to the axis of the cylinder. 

2. Heat sources and sinVa are assumed to be distributed 
on the external and internal boundaries. 

3. The temperature distribution is assumed to be deter- 
mined by the boundary temperatures and Laplace’s equation. 

4. The material behaves in an elastic manner. 

6. The variation of the elastic constants (modulus of 
elasticity, Poisson’s ratio, and coefficient of thermal expan- 
sion) with temperature may be n^lected in determining 
the thermal stresses. 

6. Plane sections initially normal to the axis of the cylin- 
der remain plane. 

7. The strains and rotation are corstant along any line 
initially parallel to the axis of the cylinder. 

The extent to which a particular structure would fulfill 
these conditions would depend on the particular circum- 
stances. The last two assumptions are appropriate when 
the cylinder is long in comparison mth its cross-sectional 
dimensions or when end conditions impose suitable restraint. 
The sixth assumption allows bending of the cylinder about 
axes perpendicular to the axial direction of the cylinder in a 
manner that might vary in the axial direction. The last 
assumption permits planes initially perpendicular to ffie 
axis of the cylinder to tahe on a warped shape as a result 
of the deformation but restricts the bending to a circular 
ai*c; that is, the radius of curvature of lines initially parallel 
to the axis of the cylinder does not vary along the length of 
these lines. The assumption of no variation of rotation in 
the axial direction is equivalent to assuming that the cylin- 
der does not twist. 

BASIC EQUATIONS 

The case of plane strain where body forces are given by 
the gradient of a potential is treated in reference 5. As 
shown by the detailed derivation in appendix A, the govern- 
ing equation of reference 6 applies to the thermal-stress 
problem defined by the preceding assumptions. This equa- 
tion is 


boundary conditions 


Two basically different types of auxiliary condition arc 
established for the set of physical conditions assumed to 
exist. These conditions are: 

(1) Conditions stemming from the nature of the forces 

applied to the surfaces 

(2) Conditions stipulating that the displacements be 

single valued 

In calculating the stresses due to the temperature distri- 
bution, other stresses such as those due to centrifugal force 
and fluid pressure are to be calculated separately and the 
total stresses are to be obtained by superposition. All in- 
ternal and external boundaries are therefore postulated to 
be free from applied forces, and the boundary conditions 
(reference 8, p. 21) become 




(4) 


where 


7 -^ 


Substitution from equations (2) into equations (4) yields 
ds"^&x5y ds 


from which 


&V dy ■ dx 
bxby ds^dx* ds 


dsby " 


ds 5x“ 


(5) 


VV=0 


CD 


along the boundaries. 



BIHAEMONIC BEIAXATION CABCTJLATION' FOB THEHMAL STRESS IN lEEEGIILAR CTLINDEHS 


73 


CONDITIONS FOR SINGLE-VALUED DISPLACEMENTS 

The defining equations for strains and rotation are 






bx 

bo 

bit , bo 


by bx 
1/bo bu\ I 

“~^bx byJJ 


( 6 ) 


From equations (6), 

bw 

bz 


= 6r 


bo 1 , 

bu 1 

bo_ 

by~^* 

btj 1 bjxt 


( 7 ) 


bz 2 bz by 

b<ri ^ ^Tty I 

by” 2 by ”^bz J 


The change in rotation for an arbitrary path P( of inte- 
gration starting at some point (a^, yo) and returning to the 
same point after enclosing (and only enclosing) the internal 
boundary Ct is 

From equations (7) and (8), the condition that <o be 
single valued is 

By use of equations (Al2) and (Al3) of appendix A, equa- 
tion (9) may be written 




bV^0 
Pi bn. 


ds= 


aE ^ bTj 

■i=i; 3 V,sr* 


( 10 ) 


The change in the z-component of displacement for an 
arbitrary path of integration starting at some point 
(^0, Vo) and retu rning to the same point after enclosing the 
internal boundary Ci is 

( 11 ) 


From equations (7) and (11), the condition for single- 
valued V. is 

£o^dy=0 (12) 

For the term involving rotation in equation (12), inte- 
gration by parts gives 

— ^ w dy=— ^£oyJ^-l-^ y dw 


where 


W“ 


0 


because of the single-valued character of to achieved by 
imposing equation (10). Furthermore, 

-1 <5 


Equation (12) can now be written 
6„dz+|7x»dy-f-|y^dz— 


In equation (13), integration of the first two terms by parts 
gives 


where 




because the strains are single valued. Equation (13) 
becomes 

(14) 

Substituting from equations (A12) and (Al3) of appendix A 
and stipulating that the stresses are single valued give 
equation (14) as 


1 — vjp, as 


b4> 


ds 


(91) 
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If the iategrations of equation (15) are performed along 
the internal boundary Ct, tlien from equations (6), 


1 

1— r 




r T- —ds^.Q 

Ci ds ox 


Because the integrand of the preceding iiitegi’al is an exact 
differential, the integral has the same value (zero) for all 
reducible paths; hence for the arbitrary circuit equation 
(15) becomes 



X 


dV^\ 
ds ) 


ds 


aE , 



(16) 


Similar reasoning from the single-valued character of the 
y-component of displacement leads to the equation 


i.(>' 




ds 


dn 


) 


ds= 


aE 


1-rJp, 


dT , dT\, 


Equations (10), (16), and (17) are similar in form to some 
equations derived by Mindlin in reference 9. The equations 
of that reference apply to integration paths taken along the 
boimdaries. 

BOUNDARY CONSTANTS 


The initial objective is to determine the stress function 
<#> so as to solve the boundary-value pz-oblem associated ivith 
equation (1). Numerical differentiation could tlien be 
performed to determine the stress components accoixling 
to equations (2). 

llie assumption of stress-free surfaces resultecl in equations 
(5). With and ae as constants of integration, these 
boimdary conditions may be integrated to give along the 
boundaries - - - 


In general, 


c)</> 






(18) 


and where aa is a constant of integration, integration of tlie 
preceding equation along a boundaiy and use of equations 
(18) give 

<f>i=anx+ai:iy+ais ( 19 ) 

along a boundary. 

Differential equation (1) and boundary conditions, expres- 
sible by equations (18) and (19), do not completely determine 
the stress distribution. No temperature terms are present, 
and the values of a«, a^, and are stUl to be specified. The 
initial step in obtaining the compatibility condition (equation 
(1)) was to raise the order of the first differential equation of 
equations (A2) by differentiating once with respect to a; and 
once with respect to y. These differentiations, although use- 
ful in simplifying the final form of the equation, require that 


additional factors be considered. In the case of multiply 
connected regions, the resulting equation (equation (1)) docs 
not preclude the occurrence of stress distributions due to 
mechanical dislocations (reference 9). Tlie occurrence of 
mechanical dislocations is eliminated by imposing conditions 
that rotation and displacements be single valued. Those 
conditions are to be used in evaluating the constants of 
equation (19). The phjzsical conditions of the problem will 
then bo satisfied, for the stress function has been defined so as 
to satisfy the equilibrium conditions; the compatibility con- 
dition is satisfied by the use of biharmonic functions; tho 
assumption of force-free boundaries is satisfied bj' conditions 
(18) and (19) on the stress function; and the constants 
a^, and are to be evaluated so as to satisfy tho conditions 
of single-valued rotation and displacements. That adjust- 
ment of the values of an, aa, and a« is sufficient to satisfy the 
conditions of single-valued rotation and displacements will 
become evident in the next section. A method of determin- 
ing the constants that will be appropriate for numerical 
techniques is needed. 

DETERMINATION OF BOUNDARY CONSTANTS BY FORMATION OF SPECIAL 

SOLUTIONS 

A method of using special solutions to determine an, o«, 
and fla for a doubly connected domain at uniform tempera- 
ture was described by Prager (reference 6). A similar 
method was suggested for the multiply connected domain 
by SouthweU (reference 10). The method of Prager is here 
extended to domains with more than one hole and with tem- 
perature distributions present. The method is then modified 
to a form that is suitable for numerical techniques. 

Let <Jiy(i==l, 2, . . . i*, . . . n; j=l, 2, 3) bo special 
solutions of equation (1) that are defined over the domain 
bounded by the external contour Co and the internal contours 
Cl, C 2 , . . . Ct, . . . On. Because equation (I) is linear, 
the products of the and the arbitrary constants Uy 
may be superimposed to give the complete solution of the 
boundary-value problem according to the scheme 

(20) 

provided that the boundary conditions for the 4>tj are properly 
selectedand that the values of the boundary constants «y are 
properly evaluated. 

The boundary conditions for the 0y must be selected so 
that thafunction given by equation (20) will satisfy equations 
(18) on all the boundaries. Because the stresses are given 
by second derivatives of <f>, the addition of a linear function 
of the coordinates to ^ will leave the stresses unaltered. Tho 
assumption is now made that this addition is accomplished 
so that the boundary constants of equations (18) and (19) 
are zero on the external boundary. The boundaiy conditions 
on the (^ij for the external contour Co are therefore taken as 

on Co 
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The boundary conditions for the on the internal bound- 
aries must now be selected so that 3n linearly independent 
solutions for the 4>n \^'iIl be obtained and so that the fimetion 
defined by equation (20) will satisfy equations (18) on all the 
internal boundaries. This selection is accomplished as 
follows: 

Let 

on all boundaries except C*. On Ct, let 


0ti=r 


dn ds 


4>ti=y 


^.= 1 ^=0 (24) 

Certain geometrical aspects of the choices expressed by 
equations (21) to (24) are now mentioned. The observation 
is made, for example, that the equation 

b4>u_ dy 
bn ds 

is a restriction on <j>ti in addition to the restriction of <f>ti=x 
on Ct, inasmuch as 4>ti=x prescribes values of only on the 
line defining the boundarj- Ct, whereas the normal derivative 

b4>ti_dy 

bn 

specifies the rate of change of as the boundary is crossed 
in a direction normal to the boundary'. 

The slope of the plane <j>a.=x in a direction normal to the 
contour Ct is b(f>tifbn; but because (t>a.=x for this plane, the 
slope may be written as bxjbn. In the direction n, the slope 
of the surface is given by equation (22) as 

b<l)ti_dy 
bn ds 


but in the ly-plane, 


dy bx 
ds~bn 


Therefore, in the direction of n, the slope of the surface ^ti 
(in equation (22)) has been taken equal to the slope of the 
plane <l>ti=x and because the intersection O't of the cylinder 
tlirough C't with the surface <pti lies in the plane ^ti=x, the 
surface <l>ti is tangent to the plane 4>ti=x along C't- The 
particular solution ati(j>tL therefore defines a surface tangent 
to the plane <f>ti=atix. Equations (23) require that <^*2 be 
tangent to the plane 4>i&=y along the intersection C"t of the 
cylinder through Ct with the surface fe, and equations (24) 


require that 4>a be tangent to the plane <^>a=l along the 
inteisection C'"t of the c^'linder through Ct with the 
surface 4>ts- 

The three special solutions, ao.4>u., and Oia^ts, asso- 

ciated with Ct are seen to be tangent to three planes, one of 
which passes through the j/-axis, one of which passes through 
the x-axis, and one of which is parallel to the xy-plane. 
Determination of a*i, Ca, and is seen to be equivalent to 
determining the slope of the plane through the y-axis, the 
slope of the plane through the x-axis, and the height of the 
plane parallel to the x^z-plane. Superposition according to 

is thus seen to satisfy the requirements of equations (18) 
and (19), but the constants Oti, att, and at* must be properly 
chosen. Superposition of aU the Oy^y will stiU leave equa- 
tions (18) and. (19) satisfied on Ct because of the requirements 
laid down on aU the 4>tj by equation (21). 

The (f>if were defined as special solutions of equation (1), 
which is equivalent to writing 

(i=l,2, . . . k, . . . 7i;J=1,2, 3) (25) 

Equations (25) together with boundary-condition equa- 
tions (21) to (24) constitute 3n boundary-value problems for 
the 3n particular solutions <^)y. These individual boundary- 
value problems with i#>y and b^nfbn specified on every bound- 
ary are now to be solved bj’ methods already described by 
Fox and Southwell (reference 7). With the assumption that 
the (j>ij are so evaluated, the next stop is to calculate the 
values of Oy. 

The method for determining the Oy may be symbolically 
expressed by substituting the complete integral as expresserl 
by equation (20) into the contour integrals of equations (10), 
(16), and (17). This substitution requires that, on each ^th 
integration path (enclosing the tth internal boundary'), con- 
tour integrals involving all the 4>u bo formed as coeflficients 
of the Cy. Formation of these contour int^rals then 
permits writing simultaneous equations for each i*th boundary 
(it=l, 2, ... n) so that 3n equations are obtained for the 
Oy where the contour integrals involving the ^y become 
coefficients of the Oy: 
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Because there are exactly 3ft equations with exactly 3ft 
unknowns (the ay), the determination of .these constants is 
sufficient to insure the occurrence of single-valued rotation 
and displacements. 

NUMERICAL TECHNIQUE 

The purposes of this section are (1) to illustrate, by a 
concrete example, the detailed steps by which a solution of 
a problem may be obtained, (2) to show, by comparison 
with a problem for which the solution is already known, that 
the application of relaxation procedures to the finite- 
difference method outlined in the preceding part of this 
report yields a method giving results that approximate the 
correct answers, and (3) to present the results of an investi- 
gation of some factors affecting the accuracy of the answers 
for some particular conditions encountered in the relaxation 
solution of a thermal stress problem. 

ILLUSTRATIVE EXAMPLE 

The detailed instructions enabling the relaxation calcula- 
tion of thermal stresses in an irregular cylinder are presented 
in appendix C. The particular problem illustrated is that of 
a concentric circular cylinder possessing an asymmetrical 
temperature distribution; however, the description of method 
may be applied to cylinders of more complex shape. The 
choice of the concentric circular cylinder enabled the relaxa- 
tion work for the stress functions to be confined to a 90° 
sector. There would be no fundamental distinction in carry- 
ing out the calculations for an irregular profile — relaxation 
would merely have to be performed over the entire cross 
section entailing more labor. (Although circular bound- 
aries were involved, the advantages of using polar coordinates 
were not utilized. The use of rectangular coordinates in the 
presence of circular boundaries involves boundary technique 
problems typical of a more irregular region.) 

The results of the relaxation calculation (reference .11) 
according to Laplace’s equation for the temperature distribu- 
tion are presented in figure 1. Results of relaxation osc- 
ulations for special solutions of the biharmonic equation are 
presented in figures 2 to 1. Contour integrals were calculated 
as illustrated by table I and the stress function is presented 
in table II. Tangential stresses (table II) were calculated 
by computing the second derivatives of the stress function 
with respect to radius (reference 8, p. 53) according to the 
five-point method of reference. 12 described in appendix D. 

Exact values of radial and tangential stress were calculated 
as indicated in appendix E and are also listed in table II. 
The maximum tangential stress is seen to be. much larger 
than the maximum radial stress. The error in the relaxation 
calculation of the maximum tangential stress was about 6 
percent. Comparison of the values of tangential stress 
calculated by the exact method with values calculated by 
the relaxation method is also presented in figure 6. 

INVESTTGAHON OF FACTORS AFFECTING ACCURACY 

Several factors influencing the accuracy obtainable in 
calculating thermal stresses by the method just described 


were investigated by applying various calculation techniques 
to a problem for which answera could be calculated by exact 
mathematical methods. The example chosen consisted of a 
concentric cylinder subjected to a symmetrical temperature 
distribution. Details of the relaxation calculation arc pre- 
sented in appendix F. 

The results of the relaxation calculation for the tempera- 
ture distribution are presented in figure 6. Because dimen- 
sions of the cylinder were exactly the same as those of the 
illustrative example previously discussed, the special solutions 
4>ii, and ^is determined in that example are usable in the 
present example. This situation illustrates an important 
feature of the use of relaxation metliods in calculating thermal 
stresses; that is, once the time consuming biharmonic relax- 
ation work for the special functions 4>a, <f>a, and 4>a has been 
completed for a given shape of body, relatively little extra 
work is required to study tlie effects on thermal stresses of 
changes in temperature distribution. 

For the symmetrical temperature distribution of figure 6, 
the boundary constants flu and were.found to vanish. 
The values of for the paths a, b, and c of figure 4 are pre- 
sented in table III. The Airy stress function and the stresses 
calculated from it by numerical differentiation according to 
appendix D are listed in table IV. 

Exact values of Airy’s stress function were calculated 
according to appendix G. The value of the arbitrary con- 
stant D in equation (G4) was adjusted to give ^=0 for 
r=12. Results are listed in table IV. Second derivatives 
of the exact values of Airy’s stress function were then cal- 
culated by, the numerical methods of appendix D to give the 
tangential stress values listed in table IV. E.xact values of 
tangentiaJ stress were calculated using the second of equa- 
tions (El) and are also fisted in table IV. 

Comparison of the errors m tangential stress of table IV 
shows that the errors associated with the. numerical differ- 
entiation of the relaxation-calculated Airy function wore 
much larger than those associated \vith the numerical dif- 
ferentiation of the exact Airy function. (Errors in the max- 
imum sti’ess were 21.5 and 6.5 percent, respectively.) This 
comparison suggests that the rela.xation-calculated Airy 
function was an important source of error. 

The relaxation-calculated stress function was calculated as 
the product aii<t>is (ajipendix F). An exact value of an 
was calculated by observing that the exact values of ^ in 
table IV range from zero at r=12 to —467,627 at r=4; 
whereas, in figure 4 the values of range from zero at 
r=12 to 1000 at r=4. The exact value of On is therefore 


-467,627 

1000 


467.627 


As indicated by the data of table HE, the errors associated 
with the processes of integrating difference quotients to cal- 
culate the boimdary constants can be significant but are 
snaall if the calculated values of the boundary constants are 
averaged over several paths. 




FiauEK 1.— Ilusulta of rclftifttlon calculation of OByzooK^trlcal tamptjraturo dlalrlbullon. Tompumturoa In ilogrina Pubronholt; on Innor bounrlaryp T— 0; on miter boundary! T— flOO-f-lOOO eos 0 , 

Tomporalure vaUioa on biU'ieraiion potha aro obtalnm) by averoiirinK Umporaturo voliiua at four naunat nodal polnla. 
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Figube 2.— Kesolts of relaxation calculation of Values of on Integration paths are obtained by averaging values at four nearest nodal points. 
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Figure 6— Resolts of relaxation calcnlatloa for temperature distribution In concentrio hollow cylinder wich symmetrical temperature distribatioa. Inner radius, 4 Inches; outer radius. 

12 inches; Inner surface, 0^ F; outer surface, S00° F. 
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The extent to which errors in the relaxation calculation of 
the Airy function influence errors in the stresses was evalu- 
ated by rounding off the previously determined exact values 
of Airy’s stress function to three significant figures. Corre- 
sponding stresses and errors are presented in table IV. In 
general, the eiTors associated with the three-figure Airy func- 
tion are seen to be significantly lower than those associated 
with the relaxation-calculated Airy function. Apparently, 
the equivalent of tlirec-figure accuracy was not achieved in 
the relaxation calculation of the Airy fimction. Compari- 
son of the values of tangential stress calculated by the exact 
method with values calculated by the relaxation method is 
also presented in figure 7. Improved accuracy could be 
accomplished by (1) further reduction of residuals with the 
introduction of another significant figure, (2) the use of a 
finer net spacing, (3) the use of more elegant finite-difference 
methods, or (4) some combination of (1), (2), and (3). A 
critical discussion of some factors influencing the accuracy 
of relaxation procedures is contained in reference 13. 


CONCLUDING REMARKS 

This investigation has yielded a numerical method for cal- 
culating thermal stress in a cooled irregular cylinder possess- 
ing one or more cooling passages under steady-state temi)er- 
ature conditions. Application of the method to structures 
such as internally cooled turbine blades is suggested. The 
use of relaxation methods and elementary methods of finite 
differences has been foimd to give approximations to correct 
values when compared with previously known solutions for 
concentric circular cylinders possessing symmetrical and 
asymmetrical temperature distributions. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, May fO, 


APPENDIX A 

derivation of biharmonic equation 


An outline of the derivation of the governing partial 
differential equation for the stress function is given in refer- 
ence 4. A detailed derivation using the conditions imposed 
in the section on assumptions follows: 

The defining equations for the normal strains (reference 8, 
p. 7) are 




«•= 


bu' 

da; 

bv 

by 

bw 

bz 


(Al) 


and the defining equations for the shearing strains are 


S'!! , da 

bu , dw 
^“=di+di- 


y|/^ 


d i) ■ dw 
’dz' by 


(A2) 


The defining equation for the rotation (reference 8, 
p. 162) is 


^1 /by du\ 
2 vda: by) 


(A3) 


The assumption that the strains and rotation are constant 
in the direction of the z-axis (axis of the cylinder) permits 
iviiting 

de^ d€„_d«,_dT^ 


bz bz bz bz bz bz bz 


= 0 


(A4) 


The conditions that the strains be compatible with dis- 
placements specified by u, v, and w are (reference 8, p. 190) 




d*6y 





dx* 

dxdy 



d**. 

b‘e,_ 

d*7„ 



dz*“^ 

dy* 

dydz 





,&*7» 



dx*"^ 

dz* 

dxdz 


d*€, _ 

_1 d 


ti , £>7 i» 

bybz~ 

~2 bx 

\ dx ‘ dy ' dz 

bhy _ 

-lA 



i_&7a*\ 

bxbz 

2 dy 

\ dx 

dy ^ 

^ bz J 

d*e. _ 

=--( 

1 

1 

* 

^ / 




(A5) 


From equations (A3) and (A4), equations (A5) become 
d*«* , d*e,_d^' 


dy* 


da:* dxdy 


dy* 

da:* 

d*t, 

bybz 

b\ 

bxbz 


=0 


=0 


= 0 


dxdy 


=0 


(AO) 
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The first of equations (A6) remains, ^rhereas the fourth, 
anil fifth of equations (A6) vanish identically. The second, 
third, and sixth of equations (A6) may be written 

^ 

dy-dz~dx^dz~dxdybz 

Because of the assumption that plane sections initially 
normal to the axis of the cylinder remain plane after the 
deformation, the displacement in the z-direction may be 
written 

v'=f(.z)x + g(z)y+Hz) 


show'ing that equations (A7) vanish identically. The pre- 
ceding proof that the last five of equations (A6) vanish 
identically is essentially a demonstration that these com- 
patibility conditions are satisfied as a consequence of the 
initial assumptions. Of equations (A6), the only nonvanish- 
ing equation is now the compatibOity condition 


5-e, b-yxy 

by- dx® bxby 


(A8) 


The generalized Hooke’s law equations (reference 8, p. 204) 
are 


1 


iz-E 


[<Xs—vi.<r,+<Tt)\-k-aT 


6,=^[<r,— v(<rx-t-(T*)] + «r >■ 
[<r»— v(o'*+<rr)] + «T 


(A9) 


Elimination of v, between these equations yields 
Cx= ^ [( 1 — 1 -|- v) V, -H 1 -1- v) a£T— piTe,] 


WAIO) 


[(1 — Hi +«')o'x+(l -{-v)aET— J 
Air 3 ’’s stress function is defined by 


■) 

^^~by^ 

b^ 

bxby^ 


( 2 ) 


Equations (AlO) and (2) are used with the direct method of 
calculating thermal stresses; that is, the equilibrium, bound- 
ary, and compatibility conditions are used without regarding 
temperature terms as body and surface forces. The equilib- 
rium equations (reference 8, p. 195) are therefore written 
without body forces. These equations were reduced to 
those of reference 8, page 21, by naing the assumption that 
plane sections initially normal to the axis of the cylinder 
remain plane after the deformation, so that 


bcr^ , 
bx^ by 




bg-f . brxit « 

by^ bx 


(All) 


Substitution of expressions (2) in equations (All) shows 
that the stress function has been defined so as to satisfy 
the equilibrium equations identically. The condition re- 
maining to be satisfied by ^ is the compatibility condition 
(AS). Substitution of expressions (2) in equations (AlO) 
yields 


E\J>y^ 


P 


P 


^^—v^4>—vEt^+(,l +y)aT 


Y (A12) 


The shearing strain is expressed (reference 8, p. 10) in 
terms of the shearing stress by 


Try 


2{l+y) 

E 

2(1 -hr) b^ 
E bxby 


(A13) 


By use of equations (A6), (A12), and (Al3), equation (A8) 
can be written 


1 — V 

For the assumed temperature conditions. 


v®r=o 


and hence the governing equation for 4> is 

V*<#.=0 


( 1 ) 


APPENDIX B 

CALCULATION OF AXIAL STRESS 


After (Tx and o-, have been calculated according to equation 
(2), the normal stress o-, in the axial direction may be deter- 
mined. The third of equations (A9) is 


from which 

o'x=exE'-i-r(vx-f 0-,) — a£’T (Bl) 

As stated in the section assumptions, the cylinder is 
free from applied forces and therefore the total force in the 
axial direction and the x- and y-components of a bending 
couple are zero. 


€i=^ k,— J'(o-x-l-<r»)]-l-ar 
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J'J' (ftdxdy =0 '| 

Jj xcr,dxdy=0 

j^^ya,dxdy=Q j 


(B2) 


Because assumptions numbers 6 and 7 in the section 
ASSUMPTIONS, the strain in the 2 -direction may be written 


e,=nx + by+c 


(B3) 


where a, b, and c are constants. 

Substitution of expressions given by equations (Bl) and 
(B3) in equations (B2) gives 

aJJ xdxdy-\-b ^^^ydxdy-\-c^^ dxdy 

[v(ffx+a’g)—aET]dxdy 


y (B4) 


ajj^x^dxdy+bjj^xydxdy+cjj^xdxdy 

= —^ ^J^x[v{ffr+ cr^)—aET\dx dy 

<i^^^xydxdy-\rb ^^^y^dxdy-\-c^^^ydxdy 

= —^^^^y[v{<Tx-^(r^—aET\dxdy 


Equations (B4) are three simultaneous equations from 
which the constants a, b, and c may be evaluated. Substitu- 
tion of the expression for given in equation (B3) for in 
equation (Bl) gives .. . . 


a,=E{ax-\-by-\-c)-\-v{<Tx+(r^ — aET 


( 3 ) 


from which the axial stress^ may be determined. 

If the origin of the coordinate axes is taken at the center 
of gravity of the cross section, 

a:da:dy=J'J^ydz(ij/=0 (B5) 


and in any case, 

^^^dxdy=A (B6) 

where A is the area of the cross section. 

Where sufficient symmetry is present, the orientation of 
the principal axes of inertia can often be determined by 
inspection. A method of determining the principal axes and 
principal moments of inertia for arbitrarily shaped areas 
is presented in reference 14. If the orientation of the co- 
ordinate axes is chosen so it coincides with that of the prin- 
cipal axes, then 


JJ^y‘dxdy=/x 

m 

J'J* xydxdy—0 

(B8> 

0* x^dxdy=Iy 

(B«> 


where 7* is the moment of inertia about the a:-axis, and 7^ is 
the moment of inertia about the y-a.xis. Simplification of 
equations (B4) using equations (BS), (B6), (B7), (B8), 
and (B9) furnishes explicit sofutions for a, b, and c as follows: 


El, ■JL ^ [v{a-x-{-(y,^ — aET]dx dy 
JJ^yl’'i<^*'^'^v)—otET]dxdy 
JX [>'(o’»+crtt) — aET]dxdy 




HBIO) 




Substitution of equations (10) in equation (3) gives 


fft=v(ifx-\-a,) — aET — 


X 

z 


IL 


(v»+0 


— aET]dxdy— 


Y JJ'^ylp(o'z'i-o-t)—ctE 2 ^dxdy— 


(Bll) 


APPENDIX C 

DETAILS OF RELA3EATION CALCULATION OF TANGENTIAL STRESS IN CONCENTRIC CYLINDER WITH 

ASYMMETRIC TEMPERATURE DISTRIBUTION 


The temperature distribution was assumed to be deter- 
mined by temperatures in degrees Fahrenheit of zero on 
the internal boundary of 4-inch radius and 600 -fl 000 cos S 
on the external boundary of 12-inch radius. The values of 
the elastic constants were assumed to be 

a=8.0X10“» (in./in./°F) 

£^=17.5X10« Ob/sqin.) 


p=0.3 

The results of the relaxation calculation for the tempera- 
ture distribution according to the technique described in 
reference 11 is presented in figure 1. In all calculations, 
the origin of coordinates was located at the center of the 
cylinder. In the present case (one internal boundary), 
the indices of equations (21), (22), (23), and (24) become. 
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i=A'=re=l 
;=1, 2, 3 

In order to eliminate the use of decimals in the work of 
relaxation, equations (21), (22), (23), and (24) were modi- 
fied as follows: 


On the exterior boundaiy, 




4>is 


^4>U 

dn 


Sn bn 


On the interior boundary. 


To that end, the appropriate contour integrals were calcu- 
lated as illustrated for functions involving i^»n and T along 
path c in table I. The numbering system for points in 
the tables corresponds to the numbering system exhibited 
by the small diagram in figure 2, where y is the station along 
the path of integration andy=l on the positive a;-axis. 

With the use of averages of contour integrals for paths 
a, h, and c, the simultaneous equations (26), (27), and (28) 
become 

1241.3 a,s=-200X2853 

8466.3 flu=0 

8466.3 «!!= -200X9532.3 


i^»ii=250 X 


^_2o0^ 

dn as 


from which 


a 13 — 459. i 


012=250 y 


bn 


250^ 

ds 


013=1000 


^0L3 

bn 


0 


Ou — 0 

aii=— 225.2 

According to equation (20), 


The boundary-value problems just defined for the bihar- 
raonic functions 0u, 0i2, and 0u were then solved by the 
tecluiiques of references 7 and 11 to yield the solutions 
presented in figures 2 to 4. (The dimensions of the cylinder 
were such that the distance between nodal points could be 
conveniently taken as 1 inch. For more general cases a 
method of handling dimensions, nodal distances, and 
derivatives, such as described in reference 7, can be followed.) 

The next step is to calculate the values of the constants 
Uii, an, and ffu according to equations (26), (27), and (28). 


0 = flll 0n + fll2 0Ii"l“®13 013 

from which 

0= — 225.2 011 — 459.1 0is 

The solution of this equation for points along the positive 
aj-axis (0=0°) is presented in table II. Tangential stresses 
were calculated from the relaxation solution for the stress 
function by taking second derivatives with respect to radius 
according to the formulas of appendix D. These results 
are also presented in table II. 


APPENDIX D 

NUMERICAL DIFFERENTIATION OF AIRY’S STRESS FUNCTION 


The purpose of this section is to present the methods used 
in approximating to the second derivatives of Airy’s stress 
function. The function to be differentiated is expressed by 

y=J{x) 

where in the present application y is Airy's stress function 
and I is radial distance on a cross section of the hollow con- 
centric cylinder. Where p is an integer, A is a uniform tabular 
interval of the independent vaiiable, and To is an arbitrary 
point from which the distance ph is measured, let 

Xp=Xo+ph (Dl) 

The values of second derivatives were calculated by the 
fivi'-point formulas of reference 12. For the various values 
of p, these formulas are: Withp=0, 

Z?Va=Y^i(35yo-104yi-|-1142/2-56y,+ llyO (D2) 
with p=l, 


D^yi=Y^i{llyo—20yi+Gyi+4:y3—y^ (D3) 

with p=2, 

yo+16i/i— 30j/2-M6yj— yd (D4) 

with p=3, 

3/0+4^1-16^2— 20i/j-l-llj/d (D5) 

with p=i, 

I?V4=j^(llyo—56i/i-f 1147/2— 104i/3-l-35yd (DC) 

Equation (D2) was used to calculate derivatives at r=4, 
equation (D3) was used for r=5, equation (D4) was used for 
r=6 to 10, equation (D5) was used for r=ll, and equation 
(D6) was used for r=12. 
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APPENDIX E 


EXACT DETERMINATION OE_JSTRESSES FOR CONCENTRIC CYLINDER WITH ASYMMETRIC 

TEMPERATURE DISTRIBUTION 


Tlie radial and tangential stresses, respectively, in the case 
of a symmetrical temperature distribution are given in refer- 
ence 8 on page 372 by 


<rr- 


a 

aETi 1 b , b^\ T 6”1 


2(i-j,)ln- 

a 


(El) 


where a is the radius of the internal boundary and b is the 
radius of the external boundary, and where Ti is the excess 
of the inner boundary temperature over the temperature of 
the external boundary. 

From reference 4, the stresses in a concentric circular 
cylinder possessing the temperature distribution 


T=A In r-bAo-bS(AnY"+C'nY“") cos 710+ 
(— -B«r“+D„r“") sinn0 
where A, B, 0, and D are constants, are 
aEA 


CE2) 


V,=' 


2 ( 1 -.) 


r, h an^-r\ 61 

In s ; In - — 

_ r 0-— a* aj 


ffg= 


aE r(6*-0(r*-a*n , r. • 

2(1 — v) [_ (0.2+ 6*)r* J ^ 

aEA 




aEAr, b.d^b^+r\ b ,"1 

2(1 — .)L”j* J*2 6*— 0* a J 

aE f-4a2r*-(3r2-&2)(r*+a2)-| . 

2(1-. )L (tt*+5V J(<?icos0+Asm0) 


The temperature distribution in °F is specified by T=0 on 
the internal boundary and by T=500+1000 cos Ion the ex- 
ternal boundary. The sti'esses will be calculated by super- 
posing stresses for a symmetrical temperature distribution 
on those calculated for a special asymmetrical distribution. 
The given temperature distribution is resolved into two com- 
ponents: a symmetrical component defined by T=0 on the 
interior boundary and 7’= 500 on the exterior boundary; and 
the special asymmetrical component defined by boundary 
temperatures of T=0 on the interior boundary and T=1000 
cos d on the exterior boundary. 

The stresses calculated for the symmetrical component of 
the temperature distribution were calculated according to 
equations (El). 


Now. equation (E2) is to be written so as to satisfy the 
boundary conditions on the asymmetrical component of the 
temperature distribution. The boundary conditions are 


at r=a 
at r—b 


T=0 

T=1000 cos e 


In order to satisfy these conditions and equation (E2), let 

n=l 

A=Ao=B^=D,t=0 
Equation (E2) is thereby reduced to 


(E4) 


cos 6 


On r— c. 


from which 


On r—b, 


from whicli 


or 


or 


from which 


T=(+,r+^‘) 

0=^+ia+^^ cos 6 





1000 cos 6 




cos B 


+.^1=1 000 


_gi+«_,000 


Oi =1000 


a* 6 


Cr- 


1000 

6*-a2 


Use of equations (E4) and of the preceding value of Ci in 
equations (E3) shows that the components of stress, due 
only to the temperature distribution specified by 7=0 on 
r=o and 7=1000 cos 8 on r=6 are 


cos 8 




1000 aE a^b (&*— r*)(r* — a*) 

2(1—.) r® b*—a* 

1000 aE a^b ia^r^—{Zr*—b^{r^+a^ 

2(1—.) r* b^—a* 


cos 8 


y (E5) 


Eadial and tangential stresses were superposed from 
stresses^alculated according to equations (El) and (E5) to 
give the exact stresses listed in table II. 
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APPENDIX F 


DETAILS OF RELAXATION CALCULATION OF TANGENTIAL STRESS IN CONCENTRIC CYLINDER WITH SYMMETRICAL 

TEMPERATURE DISTRIBUTION 


The temperature distribution was assumed to be deter- 
mined by temperatures of 0° F on the internal boundary of 
4-ineh radius and 500° F on the external boundary of 12- 
inch radius. The relaxation solution for the temperature 
distribution is given in figure 6. The values of the elastic 
constants were assumed to be the same as those used in 
appendix C, namely, 

a=8.0X10“® (in./in./°F) 

£=17.5X10* (Ib/sq in.) 

r=0.3 


fa-260!, ^=-250 g 

^> u =1000 ^=0 

These boundary values are the same as those used for the 
preceding problem involving the asymmetrical temperature 
distribution. The biharmonic functions 4>n, tpn, and <j>it are 
presented in figures 2 to 4. Calculation of the contour 
integrals by the methods described in appendix C and 
solution of equations (26) to (28) show that 


As in the asynmietrical case treated in appendix C, the 
indices of equations (21) to (24) were taken as 

i=I*=n=l 

y=l,2, 3 

and these equations were modified as follows; 

On the exterior boundary, 


<t>n = <l>ii 


^011 

£)« d/I 


^<^18 

d/l 


0 


On the interior boundaiy, 
<^>ii=250z 


^-250^ 
on as 


flu — Qii — 0 

and that the values of for paths a, b, and c are as given in 
table III. The values of Airy’s stress function as determined 
entirely by the numerical method were then calculated 
using the relaxation values of and the average of flu for 
paths fl, h, and e according to 


0 — fll8 <^13 


Values so calculated are presented in table FT. Tangential 
stresses were calculated from the relaxation solution for the 
stress function by taking second derivatives with respect to 
radius according to the fomulas of appendix D, These 
resiilts are also presented in table IV. 


APPENDIX G 


EXACT DETERMINATION OF AIRY’S STRESS FUNCTION FOR CONCENTRIC CYLINDER WITH SYM.METRICAL 

TEMPERATURE DISTRIBUTION 


The purpose of tliis section is to present the formula used 
in calculating Airy’s stress function for the concentric circular 
cylinder with a symmetrical temperature distribution. A 
possible form for the stress function is given in reference 8 
on page 55 as 

<t>=A In rd-£/-* In r-+-0*-(-2? (Gl) 


, ^(l+2lnr)-h 

2(l-.)ln|- 


ocETt 

2(1-?') In^ 



b 


fl* 



a 


(G3) 


where A, B, C, and D are constants and the corresponding 
stress components are given by 


ff,=p+£(H-2 1nr)+2C 

ffi 4+5(3 -1-2 In r)+2C 


(G2) 


Comparison of equation (G3) with tlie first of equations 
(G2) shows that necessary conditions on the constants of 
equation (Gl) are 


A= 


aETt 


a^b^ 


6 6^-fl* 

2(1 — r)Ln - 
a 



The thermal stresses are given by equations (El). For 
comparison wdth the first of equations (G2), the first of 
equations (El) is written 


B 


aETt 1 

2(l-v)ln-^ 

a 



aETt 

2(l-j>)In- 

fl 


fl*6* 


Inf 

a 


p+ 


aETt 
2(1 — r) In 



a 


a 

b-—a‘ 
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Because the stresses are calculated from derivatives of 
Airy’s stress function, the value of the constant D in equa- 
tion (Gl) may be taken arbitrarily. The values of A, B, and 
C as just determined are substituted in equation (Gl) to 
obtain 


4>~ 


aETt 


2(l-p)ln^ 


F— ^r*In jr* 

L 2 r 4 


6 *— a *\2 


6*ln r^ln^J+/>' 


(G4) 
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TABLE I— CONTOUR INTEGRALS ALONG PATH c , 


<a) Integrals Involving v**i 


s 



dt 




B 

da 

_ tvV, 

V 

■ # 

” ds 

avV/ 

1 

-118 

-118 

0 

-96 

-170 

1 

74 

6.5 

0 

407.0 

0 

0 

0 

2 

—108 

—113 

10 

-81 

—106 

24 

6,5 

66.0 

132.0 

1 

lao 

240 

S 

-106 

-108 

-8 

-88 

-142 

59 

5.6 

- 16.6 

324.5 

2 


118.0 

4 

-86 


21 

-71 

-127 

66 

5.5 

116.5 


3 

63.0 


5 

-67 

-86 

28 

-61 

-92 

41 

5.5 

164.0 

225.5 

4 

112.0 

104.0 

8 

—29 

-67 

28 

-26 

-59 

33 

6.6 

154.0 

181.6 

6 


165 0 

7 

-52 

^20 

-28 

-27 

-69 

32 

5 

- 115.0 

100.0 

5.6 

- 120.6 

1750 


-66 

■ -52 

-14 

-42 

-82 

40 

4 

- 66.0 

mo 

56 

- 77.0 


fl 

-84 

-06 

2 

-40 

-92 

43 

8 

6.0 

129.0 

6.5 

11.0 

236.5 


-38 

-84 

36 

-48 

-73 

30 

2 

52.0 

60.0 

56 

1450 

1050 

11 

—9 

-38 

29 

-18 

-19 

1 

1 

29.0 

1.0 

5.5 

159.6 

56 

12 

9 

-9 

18 

0 

0 

0 

0 

0 

0 

56 

99.0 

0 

18 

38 

9 

20 

18 

19 

-1 

-1 

- 29.0 

1.0 

6.6 

159.5 

- 5.5 

14 

64 

38 

26 

43 

73 

-30 

-2 

- 62.0 

' 00.0 

6.6 

143.0 

- 166.0 

16 

06 

64 

2 

40 

92 

-43 

-3 

-ao 

129.0 

&6 

11.0 

-235 6 

16 

62 

66 

-14 

42 

82 

-40 

-4 

66.0 

leao 

5.6 

- 77.0 


17 

20 

52 

-23 

27 

59 

-82 

-5 

115.0 

1650 

56 

- 120.6 

-1750 

18 

67 

29 

23 

26 

50 

-88 

- 6.6 

- 154.0 

181.5 

5 

14 a 0 

- 166.0 

tfl 

8 fi 

67 

28 

51 

92 

-41 

- 5.5 

-1840 

225.5 

4 

112.0 

-1640 
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21 

71 

127 

-60 

- 6.6 

- 116. 6 

308.0 

3 

63.0 

-1650 

SI 

103 


-3 

8 S 

142 

-59 

- 6.6 

156 

324.6 

2 

-50 

- 118.0 

22 

118 



81 

105 


- 6.6 

- 55.0 

1 S 2.0 

1 


- 24.0 

23 

ns 

113 


96 

170 

-74 

-5 6 

0 

407.0 

0 

0 

0 

24 

103 

113 


81 

105 

—24 

-&5 

55.0 

132.0 

-1 

150 

24.0 

28 

106 

103 

3 

83 

142 

-60 

- 5.5 

- 16.5 

824.5 

-2 

- 6.0 

118.0 

26 

85 


-21 

7 i 

127 

-66 

- 5.5 

1 X 5.5 

308.0 

-8 

63.0 

105 0 

27 

67 

85 

-28 

51 

92 

-41 

- 6.5 

154.0 

236.5 

-4 

112.0 

1640 

28 

29 

57 

-28 

26 

69 

-33 

- 6.6 

164.0 

181.6 

-6 


.1650 

2 Q 

' 62 

29 

23 

27 

59 

-82 

-6 

-115 0 

160.0 

- 6.6 

- i 2 a 5 

170.0 

30 

t )6 

52 

14 

42 

82 

-40 

-4 

- 66.0 

160.0 

- 5 L 6 

- 77.0 


31 

64 

66 

-2 

49 

92 

-43 

-8 
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_ . TABLE I— CONTOUR INTEGRALS ALONG PATH e— Coacluded 
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TABLE II— STRESS ALONG POSITIVE X-ANIS IN CON- 
CENTRIC CIRCULAR CYLINDER WITH ASYMMETRI- 
CAL TEMPERATURE DISTRIBUTION 


Radius 

fln.> 

Stream 

function 

Relaxation 
tangentfal 
stress 
(Ib/sq. In.) 

Exact stresses 
Wsa in.) 

Eiror. la 
lebxatloa 

(Ib,'»l ia.) 

Radial 

Tungeotlal 

4 

-fiw.soii 

121,012 

0 

126,988 

-5,976 

! 5 

-652,832 

65,289 

11126 

6^806 

-4i6Ii 

1 ^ 

—554,142 

26,658 

22,172 

35,409 


7 

-421561 

:^370 

20,982 

11,675 

-9,305 


-299,825 

-11,372 

17,673 

-i277 

5,095 

9 

-181» 2ffT 

-33,615 

13.496 

-30,741 

3,074 

10 

-as,27s 

-32,635 

9,013 

-3^829 

-294 

1 11 

-26,441) 

-35.852 

1,473 

-43,542 

-7.690 1 

12 

! 

0 

-33.617 

0 

—53,012 

-19,393 


■ Positive slims denote relaxation valoes larger In absolute magnitude 
than exact values: negative sisns denote reluatloa valoes smaller fa 
absolute magnitude exact values. 


TABLE III— VALUES OF a„ AND ASSOCIATED ERRORS 



au 

Error 

(percent) 


-407.627 



—555.0 

18.7 
fl X 

PArh b 

—423.8 


—416.6 

10.9 

Average of a, 6,e 

-465.1 

.5 


TABLE IV— NUMERICALLY CALCULATED STRESvSES AND ASSOCIATED ERRORS FOR CONCENTRIC CYLINDER WITH 

SYMMETRICAL TEMPERATURE DISTRIBUTION 


RatUus 

(la.) 

Exact val- 
ues of 
. Airy’s 
- str^ 
function 

Relaxation 
values of 
. Airy’s 
- stress 
fiinctiua 

Exact val- 
OPS of tan* 

gpnlhT 

stress 

(Ib/stlla.) 

Five-point method for 
tan^nUaT stre^ using 
re&xatfon values of 
Airy’s stress hinction 
Clh/sq In.) 

FlTe-iwInt metho<i for 
tangeatXal stress osin* 
ciact values of Alrr’s 
stress function 
OWsaln.) 

Flve-pufnt method for 
tangential stress using 
Airy’s function exact to 
three figures 
Cbi'sq Inj 

Stress 

Eiroc • 

Stress 

Error . 

Stress 

Error • 

4 

-467,027 

-465,100 

66.9SS 

81,357 

14,369 

62,620 

-4^388 

60.600 

-6tl00 

S 

-catsTa 

-433. 121 

36.682 

42,235 

5^704 

36^989 

4C7 

30(600 

0 

ft 

-374.551 

-367,894 

1T,2»4 

14, S« 

-2.440 

17,254 

-30 

17,600 
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i 

-291. 67S 

-281,851 

3.636 

—620 

-4,236 

^527 

-9 

3,800 

200 

8 

-»4,799 

-I95;8fl7 

-6^746 

-8,449 

1,703 

-6^766 

20 

-t;30o 

600 

V* 

-124, 507 

-118,. 135 

-15,056 

-15,969 
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-15,054 

-4 

-13,800 

-1,300 

10 

-sa. 152 

-66,277 

-21*964 

-21,588 

-37ft 

-21,965 

1 

—23; 400 

1.100 

U 

-15i<i77 

-15^813 

-27,861 

-24,845 

—3,019 

-27,883 

19 

-27,400 

-600 

12 

0 

0 

-33, 012 

-25i780 

-7,232 

-32; 776 

-236 

-26,400 

-6,600 


* Positive quantities denote approximate valoes larger in absoUite magnitude than exact values; negative quantities denote ap- 
proximate values smaller In absolute magnitude than exact values. 




